Kinetics of Particles Adsorption Processes Driven by Diffusion. 
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The kinetics of the deposition of colloidal particles onto a solid surface is analytically studied. 
We take into account both the diffusion of particles from the bulk as well as the geometrical aspects 
of the layer of adsorbed particles. We derive the first kinetic equation for the coverage of the surface 
(a generahzed Langmuir equation) whose predictions are in agreement with recent simulation results 
where diffusion of particles from the bulk is explicitly considered. 
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The adsorption of colloidal particles onto adsorbing 
surfaces is a very common phenomenon in many fields 
of Biology and Chemical-Physics. The deposition of bac- 
teria on teeth is an example of such a process. There are 
many other systems in which suspensions are in contact 
with adsorbing solid surfaces. The adsorption of latex 
spheres onto a flat solid surface is a simple experimental 
model system, on which considerable work has been done. 
An optical microscope permits the direct observation of 
structural properties of the adsorbed layer [0|2j . 

This apparently simple phenomenon involves several 
fundamental processes that govern the kinetic behavior 
of the system and the structure of the adsorbed layer. 
First, the attractive potential that characterize the ad- 
sorbing surface plays an essential role. Depending on the 
strength of this potential, adsorbed particles can either 
desorb or move laterally or remain irreversibly adsorbed. 
Second, the kinetics of the adsorption is strongly influ- 
enced by particle transport from the bulk to the surface. 
It depends, for instance, on whether gravity or diffu- 
sion dominates the dynamics of the particles in the bulk. 
Third, the adsorbing particles interact with the adsorbed 
ones in the vicinity of the surface. This interaction is re- 
pulsive for stable colloidal dispersions and, therefore, it 
is responsible for the saturation of the adsorbing surface, 
the so-called blocking effect. 

Theories have been proposed to analyze the contribu- 
tion of these effects to the kinetics of the adsorption of 
colloidal particles. A very common model is the random 
sequential adsorption (RSA) introduced to generate 
the structure of a layer of irreversibly adsorbed colloidal 
particles. In this model, the adsorbing particles are de- 
posited sequentially at random. The overlaps are for- 
biden so that, when in a trial an adsorbing particle over- 
laps an adsorbed one, this trial is rejected and one pro- 
ceeds with a new one. Clearly, RSA takes the excluded 
volume effect into account but it cannot deal with the 
transport process of the particles from the bulk to the 
surface. However, a kinetic equation for the RSA model 



has been proposed If one deflnes the coverage of a 
two-dimensional surface 9 as irR^ps, Ps being the number 
of adsorbed particles per unit area, and R the radius of 
the particles, one then writes 
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(1) 



Here Ka is a kinetic coefficient, ps is the concentration 
in the bulk and ^rsa{9) is the available area at a given 
coverage for a new particle to adsorbe. ^rsa was first 
obtained by Widom 

Some effort has also been made in constructing a model 
incorporating a description of the transport of adsorbing 
particles from the bulk. For spherical particles, heuristic 
arguments allowed Shaaf et al. jj] to describe the kinet- 
ics of irreversible adsorption for a system near surface 
saturation. When diffusion is explicitly considered, these 
authors found that the time-dependence of the surface 
coverage is changed with respect to that obtained from 
the pure RSA model (eq. (g)). Tarjus and Viot g later 
proposed an analytic (l-l-l)-dimensional model, in which 
hard disks were adsorbed onto a line, in order to give 
more insight on the previous results of Schaaf et al. From 
ref. Q , one can infer that the shape of the particles plays 
a significan role in the time-dependence of the coverage 
near saturation. Currently, for three dimensional sys- 
tems the kinetics of adsorption over all time ranges is 
only known from simulations Q 

In this letter, we propose a model taking into account 
both the diffusion of the particles from the bulk and the 
saturation of the surface due to the blocking effect. Our 
aim is primarily to derive a kinetic equation {generalized 
Langmuir equation) describing the time-evolution of the 
coverage for a three-dimensional system. Our equation 
is valid not only near saturation but also at intermediate 
times, when the surface is only partially covered. The 
results issued from our work are in good agreement with 
recent simulations |^ (Fig.l). Our secondary aim is to 
show that eq. (pi) corresponds to the adsorption of dif- 
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fusing cylinders. This confirms previous authors 'claims 
as to the crucial role played by the shape of the particles 
in the kinetics of adsorption. 

The following analysis can be applied to a large class of 
systems, provided that they meet the following require- 
ments; (i) the attractive potential, responsible for the ad- 
sorption of the particles on the surface, is large enough 
so that the adsorption can be regarded as irreversible. 
Then the particles are not observed to desorb or to dif- 
fuse laterally on the surface; (ii) the surface-particle as 
well as the particle-particle potential is of short range, 
much smaller than the size of the particles; (iii) diffu- 
sion dominates over gravitational effects in the dynamics 
of the adsorbing particles, thus the Peclet number (Pe= 
mgR/kT, where m is the mass of the particles, g is the 
acceleration of the gravity and kT is the thermal energy) 
of the system is much smaller than the unity. 

The procedure adopted here is to consider that the 
system of adsorbing particles is in local equilibrium, so 
that a local chemical potential can be defined @-||]. Here, 
the local equilibrium hypothesis is assumed since Brown- 
ian motion permits the particles to explore large regions 
in space before they get adsorbed. Thus, one expects 
that the overall adsorption process is not determined by 
the local inhomogeneities of the layer of adsorbed parti- 
cles but by its global properties, in the spirit of a mean 
field approach. Finally, for the system under scrutiny, 
the only relevant direction will be that orthogonal to the 
plane (the z-direction) since adsoption will be responsible 
for an overall density gradient depending only on z. 

Therefore, considering that the suspension of the ad- 
sorbing particles is dilute, one can write the local chem- 
ical potential as 



^(7) = kTln[p{j)] - kT\n<i>{e,j)) 



(2) 



where 7 = z/2R is the dimensionless coordinate in the 



z-direction, R being the radius of the particles. The ori- 
gin of coordinates is thus taken at a distance R from the 
adsorbing plane so that z = coincides with the centers 
of the adsorbed particles. Therefore, it is in the region 
< 7 < 1 that the interaction between the adsorbing 
and the adsorbed particles takes place. This region will 
be referred to as the interaction layer from now on. In 
eq. (^, the first term in the right hand side comes from 
the configurational contribution as in an ideal gas. The 
second term accounts for the accessible area at a given 
height 7 and at a given coverage 9, through the available 
surface function $(^,7), which comes from an average 
over the configurations of the deposited particles at a 
fixed coverage Q. Note that, for a given coverage, the 
change of the available surface function with 7 is only 
due to the shape of the particles. It is easy to see that 
the area excluded at 7 by a single adsorbed sphere is 
simply 47ri?^(l — 7^). Thus, the available surface func- 
tion can approximately be written as a function of one 
argument, according to $(^,7) — $(0(1 — 7^)). In this 
approximation an additional explicit dependence in 7 is 
ignored, due to the fact that in performing the average 
over the configurations the true radius of the particles is 
not the scaled one at a given 7. This explicit dependence 
is not too important for our purposes since the kinetics 
mainly depends on the behavior of ^(9,^) for 7 near 0, 
which corresponds to the bottle-neck for the whole pro- 
cess. We will further assume that the chemical potential 
of the bulk corresponds to that of a dilute suspension of 
particles, /is = kThi[pB]- 

Recent experimental results indicate that the ge- 
ometrical aspects of the adsorbed layer for real systems 
at small Peclet number are practically indistinguishable 
from those obtained from RSA filling rules. There- 
fore, this allows us to further assume that <I>(6',7) ~ 
^r.sa{9{1 — 7^)). Using the uniform approximation for 
the ^-dependence of ^rsa as given in ref. |0, we can 
finally write 



[1 - x{i ~ -f^)y- 



1 - 0.812x(l - 72) + 0.234x2(1 - 72)2 + 0.084a;3(l _ ^2)3 



(3) 



where x 



3, being the coverage at saturation. 



To proceed further, we will analyze the flux of parti- 
cles going to the surface. Adsorbing particles satisfy a 
continuity equation of the form 



d_ 
di 



pin.t) 



d_ 



(4) 



where we have considered that the variables depend only 
on the z-direction through 7, and on the time. Since the 
gradient of the chemical potential (0) is the thermody- 



namic force responsible for the flux of particles, one can 
write 



J(7) = -i(7) 



(5) 



where we have assumed a linear relation between the flux 
J(7) and the thermodynamic force dp,{'y)/d^, ^(7) being 
a phenomenological coefficient in the spirit of Onsager's 
approach ||l^. In this way, from eqs. (|^) and (||) we 
arrive to a Smoluchowski equation for the density of the 
adsorbing particles inside the interaction layer 
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where we have defined the diffusion coefficient D = 
kTLR^/p assumed as constant from now on. The ad- 
sorbed particles induce an entropic barrier, which de- 
pends on their shape and on the coverage, through which 
the adsorbing ones have to diffuse to reach the wall. Cru- 
cial is to note that eq. (^) differs from the Smoluchowski 
equation obtained by Shaaf et al [Q. The cross-sectional 
area in ref Q is defined only near saturation and the 
master equation used have no meaning out of this limit. 
Thus the presence of $ in the Smoluchowski equation (||) 
is not trivial. 

We will now derive the kinetic equation. We assume 
that the process is quasistationary dp/dt ~ 0. This 
means that J(7) is a constant, standing for the fact that 
no accumlation of adsorbing particles in the interaction 
layer occurs. Thus, mass conservation leads to ^(7) = J^, 
where Js is the flux of particles reaching the adsorbing 
surface. We consider here that the density of particles in 
the bulk is the control parameter and thus express Js in 
terms of ps, with the boundary conditions 



m(7 1) = MB, 
/i(7 = 0) ^ —00 



(7) 

(8) 



The first boundary condition assumes that the chemi- 
cal potential is continuous through the upper boundary 
of the interaction layer. The second stands for an irre- 
versible adsorption process. It is equivalent to imposing 
/3(7 = 0) = 0, according to eq.(^, meaning that there are 
no free particles in contact with the surface. In this way, 
we can obtain the flux of particles reaching the surface 
in terms ps by solving the differential equation 



J, = 



D 



^ ^ 1 ^' 



(9) 



with boundary conditions specified in eqs. (0) and (H). 
We then arrive at 



where 



Js 



1(0) 



D 



PbI{0) 



(10) 



(11) 



Therefore, the equation for the density of the adsorbed 
particles has the form 



dps 
dt 



-Js 



D 



PbI{0) 



This equation can be rewritten in terms of the coverage 
by multiplying both sides by ttB?' , yielding a generalized 
Langmuir equation 



- = K,pBi{e) 



(13) 



(12) 



where we have defined the kinetic coefficient Ka = Dn. 
Therefore, the main result of this letter is a equation 
describing the covering of an adsorbing surface when dif- 
fusion as well as the blocking effect are important. 

The asymptotic time-dependence near saturation is ob- 
tained by analyzing the behavior of I (9) when 9 — > Ooc- 
In this limit, this function behaves as I{9) ^ {9 — 9oc)^^^ . 
Inserting this assymptotic behavior in eq. ( p^ ) we find 
that -0OO ^ Therefore, our equation recovers 

the power law behavior near saturation B. However, 
we can go beyond the calculations of rcf. Ef because eq. 
( p^ ) gives the evolution of the coverage in all the time 
range. We have to precise at this point that at the very 
begining of the deposition process, the quasistationary 
conditions break down. Thus, our analysis is valid only 
after this initial time range which is of the order of R^/D, 
In particular, for all values of 9 /Ooo we find a remarkable 
agreement between I{9) obtained here and the results 
reported in ref. from simulations of diffusing spheres 
(Fig. 1). The underlying physics leading the kinetics of 
the system is independent of the distance h at which one 
place the particle reservoir keeping constant ps- How- 
ever, the quasistationary condition critically depend on 
this distance due to the fact that the density profile re- 
laxes with a time scale h? /D. Note, however, that in ref. 
Q the reservoir is located at a distance of the order of i?, 
being then justified the comparison of these simulations 
with our calculations. 

Let us now consider a system in which the particles 
are cylinders (radius R and heigth 2R) that diffuse with 
their axis perpendicular to the adsorbing surface. In 
this case, the available surface function $ is indepen- 
dent of 7. Thus, from eq. (p]), one straightforwardly 
gets I {9) — ^rsa{9)- Therefore, for this system we have 
found the same kinetic equation as that proposed for the 
RSA model (eq. (||)). Obviously, the assymptotic behav- 
ior for the coverage near saturation for diffusing cylinders 
gives the RSA power law decay 9 — 9ao ~ t~^/'^. Parti- 
cles of different shape can also be studied. In particular, 
one can consider revolution bodies wich offer an avail- 
able surface function of the form <I'(6'(1 — 7^")^/"), which 
would correspond to spheres for n = 1, and to cylinders 
for n — > 00. For intermediate values of n, these particles 
look like cylinders with round edges. One can then de- 
termine the kinetic equations for such objets and find the 
corresponding power law behavior near saturation. The 
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assymptotic behavior depends on n and has the form 

e-Ooo- i2n/(l-4«)_ 

Our analysis, therefore, gives an interpretation of the 
previous simulation results pointing out the collec- 
tive effect of diffusion, blocking effect and the shape of 
the particles, on the kinetics of adsorption. Due to the 
generality of the basic statements underlying our theory, 
one can envisage its use to the study of the kinetics of 
more complex systems as, for instance, suspensions of 
elongated particles, provided that the structural proper- 
ties of the adsorbed layer are given. 
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